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Abstract 

A sum rule due to Das et al. is reanalyzed using a euclidian space approach and a 
Pade resummation procedure. It is shown that the result is essentially determined by the 
matrix elements of dimension six and dimension eight operators which have recently been 
measured by the ALEPH collaboration. The result is further improved by using the vector 
spectral function which must be extrapolated to the chiral limit. This extrapolation is 
shown to be reliably performed under the constraint of a set of sum rules. The sum rule 
is employed not as an approximation to M w + — M n o but as an exact result for a chiral 
low-energy parameter. A sufficiently precise evaluation provides also an estimate for a 
combination of subleading electromagnetic low-energy parameters. 



00 



43 

o 1. Introduction 

<D '. 

Chiral perturbation theory (e.g. [Q for a comprehensive review) is now claiming to 
reach such a high degree of accuracy in some situations that it is becoming necessary 
to deal quantitatively with radiative corrections in low energy processes. An important 
example is the pion-pion scattering amplitude for which the two-loop contribution has 
recently been evaluated [0]|3|. The relevance of this reaction for probing experimentally a 
basic issue in the spontaneous breaking of chiral symmetery in QCD is discussed in some 
detail in ref.||. Calculations of radiative corrections have started to be performed both for 
the pionium atomic bound state (e.g. [Q] and references therein), in view of an experiment 
planning to form pionium atoms at CERN[||, and for the scattering amplitudeQ[Q|. The 
framework for performing such calculations is a natural extension of the conventional chiral 
expansion to include the photon as a dynamical quantum field ||. This extension brings 
in a set of new, a priori unknown, low-energy constants. At chiral order two, a single 
constant appears (which will be denoted by C below) while at the next chiral order, one 
has to deal with fourteen new constants called k{ in the case of the SU (2) x SU (2) chiral 
group 00. 



The purpose of this paper is to reanalyze the classic sum rule of Das et al.Q. Since 
experimental data has started to become available from r decays into hadrons, the sum 
rule was discussed several times in the literature [13| []14|] [15| using as input experimentally 



measured vector and axial-vector spectral functions. One must be cautious, however, that 
the sum rule can at best provide an approximation to the 7r + — 7r° mass difference if 
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the spectral functions are not extrapolated to the chiral limit. Indeed, the derivation is 
made in the limit, m u = rrid = and, stricly speaking, the integral diverges if one uses 
physical spectral functions over an infinite range. In modern context, the sum rule must 
be interpreted as an exact result for the low-energy constant C. This constant appears in 
the leading term of the chiral expansion of the ir + — tt° mass difference 

Ml + - M 2 = ^f+ O(e 2 M 2 w0 ) + 0((m u - m d ) 2 ) + 0(e 4 ) . (1) 

The corrective terms O(e 2 M 2 ) and 0((m u — md) 2 ) involve a number of low energy con- 
stants ki and one 0(p A ) constant (Z7) respectivelyR[10]. The order of magnitude of 



low-energy constants such as ki is known from rather general considerations on effec- 
tive theories]ll[ to be ki — F 2 /A 2 , where A is the typical mass of the massive states, 
not included in the effective theory, i.e. A ~ M p (or M K , M v in the 577(2) x SU{2) 
expansion). This enables one to estimate that the corrective terms in eq.(|l]) could be as 
large as 20 — 30%. Our claim is that by a clever use of r-decay data recently released by 
the ALEPH collaboration [12] JO]] it is actually possible to perform the sum rule evalua- 
tion of C in such a precise way as to actually provide an estimate for the combination of 
low-energy constants involved in the corrective terms in eq. ([!]). 

In practice, we advocate an approach in which one first constructs the QCD correlation 
function < VV — AA > in the chiral limit in euclidian space, an idea which was proposed 



in ref . [ 16 . A key ingredient for this construction is the experimental measurement by the 
ALEPH collaboration |l5| of the vacuum matrix elements of the dimension six and dimen- 
sion eight combination of operators which control the first two terms in the asymptotic 
expansion of the chiral correlator. In euclidian space, far from the resonance region, this 
asymptotic expansion is expected to be accurate down to rather low momenta values, say 
p ~ 2 GeV. The task is then to interpolate a smooth function of p, the value of which is 
known at zero (in terms of F n in the chiral limit), in a finite momentum range. In this 
approach, the momentum integral in the [0, 00] range can then be performed exactly. It 
will be argued that the only knowledge of the two operator matrix elements (together with 
F n ) constrains the value of C to a level close to 10%. The estimate will then be refined by 
using more detailed experimental information on the vector and the axial-vector spectral 
functions. 

2. Description of the method 

The starting point is the sum rule derived by Das et al.|| (a quick derivation can be 
found in ref.[[l7j) written as an integral in four dimensional euclidian space. Performing 
the angular integration, one expresses the constant C as a one dimensional integral 



C = — 7T [ ds s 

4 16^ 2 J 



n A (s)- n v (s) (2) 



where 1X4 and fly are defined as the limit when m u = = 0, e 2 = of the form-factors 
11^4 and ily associated with the axial-vector and the vector two-point correlation function. 
Ily, for instance, is defined as 

i [ d 4 xe l P x < 0|7%(x)Vj(0)|0 >= (p^-pV) n v(p 2 ) +pV n °vV)> (3) 
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being the charged vector current V^(x) = u(x) r y fl d(x). An exactly analogous definition 
holds for IXi. 

Formula (||) is exact provided chiral symmetry is spontaneously broken in QCD with 
two massless quarks. It is of interest to further consider the SU(3) x SU(3) chiral limit 
obtained by sending m s to zero as well. However, as will be seen in the sequel, the uncer- 
tainties involved in this extrapolation are too large and do not permit a useful evaluation 
of Cq = lim ms= o C. Convergence of the integral in (|Tp follows from applying the operator- 
product expansion[18|[19|. The operators must belong to the (3,3) representation of the 
SU(2) x SU(2) group. In the limit m u = = e 2 = the only such operators that one 
can construct are of dimension six or more. The spontaneous breaking of chiral symmetry 
then implies that the vacuum expectation values must be non- vanishing. The following 
asymptotic expansion thererefore holds, 



lim U A (-p 2 )-U x 



-P 



-^6 



+ ... 



p u 



(4) 



In QCD, A6 and Ag are not exactly constants except at leading order in a s . At higher 
orders, corrections carry logarithmic-type p 2 dependences. In the leading-logarithmic ap- 
proximation, this p 2 variation is found to be rather slow, such that the approximation of 
taking constant values for Xq and As will be accurate in a reasonably large energy region. 
We will return to this point in sec. 4. The parameters A6 and As have been determined 



experimentally using r decay data|15]. The method consists in using the analyticity prop- 
erties of the two-point functions together with Cauchy theorem, which leads to equations 
like 



Jo 



dsP kl (s)[p A (s)-p v < 



dzP k \z) [U A (z)-U v (z)} 



--Ml 



where P kl (s) can be any polynomial. A convenient set is [ 20 1 



(5) 



P kl (s) 



1 



s 



2s \ 

1 + 



k,l > . 



(6) 



For k = I = the left hand side of eq.(||) reduces to a difference of total r decay rates. 
These polynomials have the further merit to suppress the contributions which are close 
to the cut in the integral over the circle so that one can use asymptotic QCD expansions 
with some confidence in the righthand side of eq.(||). Using this method, the ALEPH 
collaboration has determined the value of the following dimensionless integrals involving 
Ag and As, 



£(2n) 



4-7TZ 



M 2 J\ Z \ =M ? 

-0.058 ± 0.006 and 5^ 



dzP 



00/ 



A2n(z) 



n 



3, 4 



(7) 



to begD, ^ (6) = -0-058 ± 0.006 and 5^ = 0.0170 ± 0.0014. Ignoring the p 2 dependence 
of Ag and As (the validity of this approximation in the present situation can be checked 
explicitly for A6 and will be found in sec. 4 to be excellent) one deduces, 



A fi 



(7.58 i 0.80)10 -3 GeV D A 



-1.07±0.12)10~ 2 GeV 



(8) 
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The result for Xq is in reasonable agreement with that obtained earlier pi]]. An analysis 



of the ALEPH data, making use of negative moments performed very recently[22] leads to 
values compatible with (^) although slightly smaller. It is perhaps important to stress that, 
even though A6 and As control the expansion of a chiral limit correlator they are effectively 
correctly determined from data in which m u , ^ 0. This is because quark mass effects 
are properly taken into account in the fit as they occur in the operator-product expansion 
via operators of lower dimensionality and the contribution of dimension six linear in the 
quark mass ( involving the so-called mixed condensate) happens to vanish at leading order 
in a s []23|]. In other terms, the chiral correction to Xq is strongly suppressed. A priori, there 
is no reason for a similar suppression to hold for As, but this parameter is of lesser practical 
importance in the calculation. A key assumption which is made in the above determination 
of A6 and Xs concerns, of course, the validity of truncating the asymptotic expansion (Q) 
at order eight for p = M T . We will see below that this assumption is internally consistent 
but it is not easy to estimate the error induced by this truncation. For this purpose, 
one should be able determine more asymptotic parameters and check the stability of the 
determination. 

Let us now explain the method for evaluating the integral in eq.(||). We first split the 
integrand in two parts 

u A (-,)- ii v (s) = u e ^ v (-s) + njTvHO > ( 9 ) 

where TLj£_y(—s) is constructed from an experimentally measured part of the vector and 
the axial-vector spectral functions. We will proceed in three successive steps of approxi- 
mation, including more and more experimental information in this part, and then check 
the stability of the result. In the first approximation, we include solely the pion pole part, 

2F 2 

U e ^ v (—s) = (1 st approximation) (10) 

where F is the pion decay constant F n — 92.4 (MeV) extrapolated to the chiral limit. The 
remainder part in eq.(^), IF^^y, is reconstructed from its asymptotic expansion assuming 
that four terms in this expansion are known 

umjrr v (-s) = ^ + ^ + ^ + ^ + ... . (11) 

s^oo s S Z S d S 4 

For instance, in the first order approximation corresponding to (fiCfr), one would have 
ci2 = — 2F 2 , a± = 0, a$ = Xq, as = Ag. The point is that, firstly, we expect this asymptotic 
expansion to become numerically accurate at rather low values of the momenta, \fs ~ 2 
GeV. Secondly, the function W^ v (—s) is expected to be a perfectly smooth function 
down to s = 0. In the first order approximation, it has a logarithmic chiral singularity at 
s = with a small numerical coefficient, 

I 

limn^ e !!V(- s ) = 24 2 logs + cstt (I s * approximation) . (12) 

In higher order approximations this singularity will be exactly included in H e ^ v (—s) 
and the remainder part will be finite at s = 0. It is plausible that a simple rational 
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approximation should be able to interpolate rather precisely the remainder function in 
the range y/s = [0, 2] GeV. Imposing finiteness at s = and matching to the asymptotic 
expansion selects a unique kind of Pade approximant, 

(13) 

The parameters of the approximant being related to those occuring in the asymptotic 
expansion eq.(|i~T|) by the simple relations 

2 

a 6 — a^ag a 2 a 8 — a 4 a 6 , 
a = -5 , c = — 5 , = 04 + 020, a = a2 . (14) 

a| — 0206 «| — ^2^6 

In the second level of approximation we include into H e ^ v (—s) the most significant part 
of the 2ir spectral function together with the one-pion pole which was considered before 

oz?2 ,m? ° t T \ 

n?!v(-s) = — - dx (2 nd approximation). (15) 

s Jo x + s 

In this approximation, the logarithmic singularity (|l2"|) is properly taken into account 



provided the spectral function is correctly normalized at the origin: p 2v (0) = l/247r 2 @. 
The construction of the chiral limit spectral function p 2n knowing the experimentally 
measured one P2-k is not a completely trivial matter and will be explained in the next 
section. The remainder piece is constructed as a Pade approximant as before except that 
the asymptotic expansion parameters a» which enter are now given by 

a 2 = -2F 2 + I , a 4 = -Ji, a 6 = A 6 + /2, a 8 = X 8 - h (16) 

with 2 

I n = dxx no p 2 , (x) ■ (17) 





One can of course think of continuing in this way and include more and more exper- 
imental information such that the remainder function will become numerically smaller 
together with the uncertainty associated with the Pade interpolation procedure. The next 
step, then, would be to include explicitly the contribution from the three pion component 
of the spectral function, 

n e *v v {-s) = ^ - ^f^-dx + ^f^-dx (3 rd approximation). (18) 

What prevents one from pursuing this construction further lies in the difficulty of perform- 
ing the chiral extrapolation, which increases with the pion multiplicity. It will fortunately 
appear that convergence is very fast, such that one hardly needs to go beyond the second 
approximation . 

3. Chiral limit extrapolations 
3.1 F w 
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Extrapolation to the chiral limit of F n can be performed fairly easily using known 
results from chiral perturbation theory. The value of F, corresponding to m u = = 0, 
m s 7^ is related to F w at one loop order by the following expression 1 10 1 

^t-i^f -f^K + cW) , (19, 

which involves the scalar radius of the pion. Using for this quantity the updated value as 
given in ref.|24]: (r 2 )g = 0.60 ± 0.05 fm 2 , one obtains 

F = 86.7 ± 0.6 ± 0.5 MeV (20) 

using F % = 92 A ± 0.3 MeV[25|. The second error in the value of F is a naive order of 
magnitude estimate of the size of the O(M^) correction in eq.(|l~9|). This relatively precise 
extrapolation is to be contrasted with the situation in which one would be willing to 
further extrapolate to m s = 0. Let Fq be the corresponding limiting value of F n , it is 
related to F by the following relation pGfl 

F = F(l-^Ll(rn s B) + O(m 2 s ) ] ) , (21) 

where B is proportional to the quark condensate in the chiral limit, B = — < uu > /F 2 . 
This relation involves the low-energy constant L4. Unfortunately, there is no independent 
way of determining L4, which appears here multiplied by a large numerical factor. 

3.2 p 27r 

Let us now discuss the two-pion component of the vector spectral function. The 
possibility of performing a reliable extrapolation here is tied to the time old observation 
of vector meson dominance of the pion electromagnetic form- factor Fy. Defining p2 n in 
terms of Fy, 

3 

1 / s — AM 2 \ 1 

p27r{s) = e(s-4M 2 )^(-—^\ \F v (s)\ 2 , (22) 

an excellent fit to the data can be performed up to the tau meson mass, with a Breit- 
Wigner function for the p resonance and only a small admixture of higher mass resonances, 



with 



M 2 . M 2 ( s - AM 2 \ 3/2 



This type of parametrization guarantees that Fy(Q) = 1 and was proposed in ref.[27|. We 
will be using the numerical values obtained from a combined fit of the ALEPH r — ► 2tt 
decay data and the e + e~ — > 7r + 7r~ data[12] 

M p = 773.4 ± 0.9, Y p = 147.7 ± 1.6, p = -0.229 ± 0.020 (25) 
M p , = 1465 ± 22, T p , = 696 ±47, 7 = 0.075 ± 0.022 
M p n = 1760 ± 31, r p ,/ = 215±86. 
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Other variants in the functional form of the Breit-Wigner function B p (s) may be used 
which would result in somewhat different values of the parameters (|25|). In particular, the 
form due to Gounaris and Sakurai|^8| has better analytical properties and can approx- 
imately correctly reproduce the cut of Fy(s) in the chiral limit while the simpler form 
( B4j ) produces no cut at all. Nevertheless, for the problem at hand, we found numerically 
insignificant differences in using either parametrization. 

It is clear that extrapolation to the chiral limit will dominantly affect the lower energy 
part of the spectral function. Furthermore, the uncertainties in the parameters of the 
higher resonances p' , p" are larger than the effect of setting m u = = 0. Therefore, in 
order to obtain the spectral function in the chiral limit it is only necessary to evaluate the 

o o 

extrapolated values of the p-meson mass and width, M p and Y p . Let us now discuss this 
issue. 

In the case of the mass, firstly, one can perform a chiral expansion. At leading order, 
linear in the quark masses, the p and K* masses are expressed in terms of two independent 

o 

parameters (besides M p ) B\ and B$, 

M p =M p +2mB$ + 2mBi M K * =M p +(m s + fh)B s + 2mB 1 . (26) 

o 

One needs in principle to know both of these parameters in order to deduce M p . The 
parameter B% is easily obtain K* — p mass difference, 

2(M K * - M ) m s , . 

2mB 8 = ^-^ p —, r = -^~ 26 (27) 

r — 1 m 

(using the standard chiral expansion framework for evaluating the ratio m s /m). The value 
of Bi, on the other hand, cannot be simply determined, but this parameter is suppressed in 

o 

the large N c limit and thus should be smaller than Bg. Neglecting B\ gives M p —M p ~ —10 
MeV. The chiral expansion of the vector meson masses has been pursued recently beyond 
linear order [29] [30]. Including the leading correction, which are of order 0(mg /2 )|| gives 



for p-meson mass in the form, 



o q 2 M 2 o 4 a 2 M 3 

M P = M f + 2MBs + Bl ) - fL-j| M K p + ) - JL^ . (28) 



Here, the parameter g can be determined approximately to be g ~ 0.60 [31] [30] and 

o 

Mk= lim mii ,m d =o Mk- The parameter B% can, again, be determined from the K* — p 
mass difference and one finds that its numerical value is essentially the same as in the 
linear expansion. The corrective terms, even though suppressed in the large iV c limit, turn 
out to be relatively large and approximately cancel the contribution proportional to B%. 
Further corrections of order 0(m 2 ) are also generated at one- loop which were computed in 
ref.|3(|. This contribution depends on a rather large number of parameters. We will not 
attempt to take it into account quantitatively but simply use the qualitative fact that it 
goes in the sense of reducing somewhat the large effect of the 0{m q ) contribution, such 
that one can estimate with some confidence that the p mass in the chiral limit should lie 
in a range, 

-10MeV<M p -M p <0 . (29) 
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Concerning the chiral limit of the width of the /9-meson, we may also try to follow a 
similar approach and expand to linear order in the quark masses. Unfortunately, even at 
such a low order and dropping Zweig rule violating terms, there still remains too many 
undetermined constants. The most general chiral lagrangian terms describing vector meson 
coupling to pseudo-Goldstone boson pairs (using notations as in ref.|33|]) linear in the quark 
mass matrix are 

C vpp = % (< > +71 < {x (+) , > +72 < V^x i+) u u >) . (30) 

Note that wave-function renormalization effects of either the chiral fields or the vector 
meson fields can effectively be absorbed into the parameter 71. It turns out not to be 
possible to determine the three constants Gy (which determines the chiral limit width) 
and 71, 72 independently. Qualitatively, at least, this approach suggests, from the phase- 
space factor and the pion momentum dependence of the decay matrix element, that one 
should expect an increase of the p-meson width in the chiral limit of the order of 20%. 
This is a rather large effect and it must be properly taken into account. 

As a way out of these difficulties, one may construct a set of sum rules involving the 
difference of the spectral functions py — py. To the extent that the lower part of the 
integration region dominates, such sum rules will efficiently constrain the chiral limit of 
the p- meson parameters. One derives a first sum rule by considering the combination of 
Uy(—s) minus its chiral limit counterpart ITy (— s). Asymptotically, one has (e.g. pjj), 



lim s (Uy(-s)- U v (sj) = (31) 

s — ►OO \ / 

lim + ^M)( mu ( s ) + m d {s)f + (1 + ^M)(m u ( s ) - m d (s))< 

s-^co Sir z I IT TV 

Hence, using a spectral representation, the following sum rule must hold, 

dx (py(x)- f>y (a?)) = . (32) 



A second second sum rule, with even better convergence properties, is obtained by con- 
sidering the following s = limit, 

n v (o)- te (flv(-) + ^iog^) (33) 

This expression can be evaluated in two different ways. Firstly, one can use the chiral 
expansion of the vector correlation function: a very good level of precision can be reached 
thanks to the calculation at two-loop order by Golowich and Kambor[|35|]. Secondly, one 
can write down a spectral representation: here it is convenient to split the integration 
range into [0, 4M n 2 ] and [4M„- 2 ,oo]. In the first range, the integral can be performed 
explicitly, using the one-loop expression for the spectral function p v . Equating these two 
evaluations, one derives the second sum rule, 

Pv{ X )-Pv 0?0 = 1 f 2 A 

4M2 x 127H V 3/ 
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Ml - 8 3 .- - . r p 2 1 p 2 1 

+ 288^ l v /6 - l0g4+ 3 + 2 (Z5 - Z6 H 1 °^ + 4 l0g M|-4 

8 ^[Q(^) + 2^) + ^^] + 0(A^). 



F 2 



In this expression, and Iq are low-energy constants which appear at 0(p 4 ) []!(]] and which 
are well determined, while R(p 2 ) and Q(p 2 ) are 0(p 6 ) constants p5| (the appearance of 
Mk in the above expression is related to the fact that these constants are appropriate 
for the three- flavour chiral expansion). One expects R(p 2 ) to be suppressed compared to 
Q(p 2 ) because of the Zweig rule (for values of the scale p of the order of 1 GeV) and the 
latter constant was evaluated from a sum rulc|3f 



Q(M 2 ) = (3.7 ± 2.0)1(T 5 . (35) 

This enables one to evaluate the entire 0(M 2 ) contribution on the right-hand side of 
eq.((3j). The set of two sum rules ( |32"1) and @) can be considered as a set of non linear 

o o 

equations from which one can determine M p and T p . We have analyzed this system 
numerically, and found that it has a solution, which is unique in a physically meaningful 
range. Corresponding to the central values of the parameters cited above and including 
only the two-pion component of the vector spectral functions, one obtains 

M p -M p = -2.4 MeV T p = 180.8 MeV . (36) 

The uncertainties in this result come from two sources. Firstly, there is an uncertainty 
in the integrals of py coming from experimental errors in the parameters describing py. 
Varying the parameters in ( ^5|) one finds that, essentially, the error on M p and T p are the 

o o 

only ones that matter and that they translate into identical errors on M p and T p . Secondly, 
we have neglected in the integrals the contribution of components in the vector spectral 
function other than 2ir, i.e. py ', py K , p*y , ■ ■ . Evidently, one expects the first of the sum 
rules to be more sensitive to these contributions which set up at higher energies. One 
can make a rough estimate of the influence of these components using the quark-hadron 
duality idea, i.e. modelling the sum of all contributions by a continuum, 

P T\s) = ^e{s-M 2 ont ) (37) 

normalized to the asymptotic QCD prediction and starting at some threshold mass M con t- 
A typical value used in sum rules analysis is M cont ~ 1.5 GeV. For the problem at hand, 
we need to know also how this continuum mass varies when going to the chiral limit. 
There is of course no way to precisely evaluate that, but it seems not unreasonable to 

o o 

assume — 10 ^ M cont — M cont ^ 10 MeV, which leads to a variation A M p = ±8 MeV. The 
conclusion is that the chiral mass is, in fact, not determined to a better accuracy from 
the sum rules than it was from the chiral expansion as discussed above. Imposing that 

o 

the sum rule result be the same as that found before, i.e. M p —M p = —5 ± 5 MeV, is 

o 

achieved by taking the continuum mass parameter in the range M con t —M con t = 2 ± 6 



9 



MeV. Solving the two sum rule equations simultaneously yields the chiral mass and width 

o 

as approximately linear functions of M cont —M cont and they are found to lie in the range 
M p —M p = — 5 ± 5 MeV Y p = 180.0 ± 1.5 MeV . (38) 

One observes that the width gets determined with a much smaller error than the mass. 
The spectral function p2 n and its chiral extrapolation are shown in Fig. 1. From this figure 
one observes, in particular, that the influence of setting m u , rrid to zero is felt mostly in 
the low-energy region, y/s < 1 GeV, consistently with the starting point assumption. 

0.14 



0.12 
0.1 
0.08 

CM 

0.06 
0.04 
0.02 





0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 

s 1/2 (GeV) 

Figure 1: Two-pion component of the vector spectral function and its chiral limit extrapo- 
lation obtained from solving the non-linear system of two sum rule equations, as explained 
in the text. 

3.3 p 3w 

The spectral function piece p^-w is not known to the same accuracy as p2w Furthermore, 
it will appear that extrapolation to m u = = is plagued with larger uncertainties. 
However, p 3n is not an essential ingredient, its explicit inclusion turns out to have very 
little effect and only serves to verify the stability of the calculation. For this purpose, 
an approximate knowledge of p 3w may be sufficient. As before, one expects a sizable 
contribution from a resonance, the ai(1260) in this case. However, because the a\ has 
a larger mass than the p and especially because it has a much larger width it is more 
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questionable that the background contribution will be negligible. We will anyway follow 
the model of Kiihn and Santamaria f27fl which assumes complete dominance of the a\ and 
matches with the correct chiral 0(p 2 ) behaviour of the axial current matrix element at 
low energy ( note that the 0(p 4 ) expression has been recently worked out |37|1 ). One 
assumption in this model is that the a± decays via a two step process: ai — > pir — > 3ir 
or ai — > p'n — ► 3-7T with a small probability . In principle, nothing prevents the a\ decay 
to proceed also via the a\ — * air channel^. A clear signature for this process would be 
a difference in the a\ decay rates into 2ir~it + and 27r°7r _ . These two rates have now 
been measured separately for the first time by the ALEPH collaboration fi"5j| and found 
to be equal to a very good precision (R—+ = 9.1 ± 0.2%, R°°- = 9.2 ± 0.2%). This 



measurement supports the decay model of ref . 1 27 ] - This model is embodied in the following 
parametrization 

< 7r-(q 1 )7r-(q 2 )7r + (q 3 )\u^ 1 5 d\0 >= -i^B ai (s) (B pp ,(s 2 )Vf + B pp ,{ Sl )V£) (39) 



with 



and 



Vr = q?-q$-Q» Q - {<li Q = qi + q2 + q3, « = (Q - qif 



B p ( Si ) + /3 f B p ,( Si ) M t 



(40) 



BpASl) ~ 1 + ' BaAs) - Ml-s- l M ai T ai g(s)/ 9 (Ml) (41) 

where g(s) is a three-body phase-space integral which must be computed numerically (see 
ref.|27j for more details 0). We have determined the a\ mass and width as well as the 
decay parameter j3 from a simple-minded fit of the ALEPH data^ assuming energy 
independent errors. The resulting values for the a\ parameters obtained in this way are, 

M ai = 1.28 ± 0.01 GeV, T ai = 0.67 ± 0.05 GeV, = -0.27 ± 0.03 . (42) 

The experimental invariant mass distribution for the mode r — > ir~ tt~ ir + v is shown in 
Fig. 2 together with the result of the fit using the above parametrization. 

Now we would like to construct the chiral limit extrapolation of the 3n spectral func- 
tion. As before, we disregard the modification of the parameters associated with the p' as 
it makes a relatively minor contribution to the spectral function. Concerning the p meson, 
the extrapolation of its mass and width were discussed in the previous subsection, there 
essentially remains to estimate the modification of the a\ mass and width parameters. 
Concerning the mass, one encounters the first difficulty that the quark mass matrix not 
only shifts the 1 ++ multiplet but also mixes the states with non-zero strangeness with 

1 We approximate, as usual, a strongly interacting pion pair in an S-wave by a fictitious or real but very 
wide a meson. 

2 An approximate analytical form of g(s) is given in this reference but one must be careful that it is 

only valid for the physical value of M n an d becomes incorrect for Mg = 0. 

3 The data can be found on the website http://alophwww.cem.ch/ALPUB/paper/paper98/! . 



11 



0.006 




s (GeVl 



Figure 2: Branching fraction for the mode r — > ■k~it~it + i> as a function of the three-pion 
invariant mass squared. The experimental results from ALEPH are displayed together 
with our ht based on the Kiihn-Santamaria parametrization. 



those of the 1 + multiplet. Expanding to linear order in the quark masses, assuming ideal 
mixing, and using the ss member of the multiplet gives, 

M B ~M ai - M ^ 510 \- Mai ■ (43) 
(r - 1) 

This estimate must be considered as very approximate because of the additional problem 
that the assignment of the /i(1510) as the ss member of the a± nonet[^] is far from 
certain |38|1 . 

The value of the a\ width, finally, in the chiral limit is constrained by a sum rule 
exactly analogous to ( p2| ) in the axial channel, 

roo 

/ dx (pa(x)- p A (x)) = 2F 2 - 2F W 2 . (44) 
Jo 

In this equation the one-pion component is excluded from pa and its contribution appears 
on the right-hand side. As before, the additional assumption must be made that this 
equation constrains mostly the low energy part of the spectral function and, as a conse- 
quence, can essentially be interpreted as an equation for the 3tt component of pa- Using 
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M ai -M ai = -10 MeV, the sum rule gives the chiral width: T ai = 0.70 GeV. We shall 
be content with a single sum rule here even though it is possible in principle to exploit 
a second sum rule in analogy with the case of the p meson. The result for the physical 
three-pion spectral function and its chiral limit is displayed in Fig.3. 

0.06 



0.05 



0.04 



K 0.03 
cn 
Q. 

0.02 



0.01 





0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 

s 1/2 (GeV) 

Figure 3: Three-pion component of the axial-vector spectral function and its chiral limit 
extrapolation. 

4. Results 

Now that we have defined an approximation scheme for calculating n^_y and EP^^y 
it is straightformward to compute the sum rule integral, eq.(||). Before we do so, it is 
instructive to have a look at the integrand, which is displayed in Fig. 4 and Fig. 5 for the 
three levels of approximation. 

Fig. 4 shows the low energy region < y/s < 2 GeV. One might believe that this part 
will dominate the integral, it actually turns out that the asymptotic tail makes a non 
negligible contribution of approximately 20%. It is one advantage of this method that it 
introduces no error due to truncation of the integral. One observes that approximations 
2 and 3 generate curves which can hardly be distinguished. Fig. 5 shows a region of larger 
values of the integration variable s from which one can appreciate the approach to the 
asymptotic regime. The two-terms asymptotic expansion is seen to be accurate at the 10% 
level for \fs = M T and becomes very accurate provided \fs > 2.5 GeV. While the sum 
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Figure 4: Integrand of the Das et al. sum rule in the three successive approximations. 



of Ifjj^y and li^-v appear to be remarkably stable they are individually quite different 
from one approximation to the other. This is illustrated in Fig. 6 showing U r A e ™y, which 
is the part where the Pade interpolation procedure is used: the figure shows how this 
part becomes smaller as one includes more experimental information from the spectral 
functions. The curves are seen to be smooth, flat, an exhibit no change of sign thereby 
justifying, a posteriori, the use of a simple rational approximation. 

We can now perform a stringent test of the interpolation procedure by considering the 
integrand at low energy, comparing it with the chiral perturbation theory expectation, 

= h ■ (45) 

The low-energy constant Z5 is known from the one-loop analysis of the charge radius of 
the pion (r 2 ) v and of the pion radiative decay amplitude, it — > ewy: I5 = 13.1 ± 1.3 (using 
F w = 92.4). An additional 5 — 10% uncertainty is expected from 0(M%) contributions 
to these observables, which necessitate a two- loop analysis (at present only ir — ► ez^7 has 
been analyzed at this level of accuracy^]). Using our construction for Ua-v(— s) and 
computing numerically the limit (fl5|) we find, 

k = 11.8 (46) 
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Figure 5: Same as Fig. 4. Also shown for comparison is the two-terms asymptotic expansion 
of the integrand. 



(the result differ in approx.2 and approx.3 by less than 1%) which is slightly smaller but 
compatible with the one-loop determination quoted above. This is a rather non trivial 
check of the quality of the interpolation from the low-energy domain of the chiral expansion 
up to the domain of large energies, where the operator-product expansion makes sense. 
The result of the sum rule evaluation of the low-energy constant C are displayed in Table 
1 below: we show firstly the dimensionless quantity Z = C/F 4 which is of order unity and 
also the quantity 5 = (M 2 + — M 2 — 2e 2 C / F 2 ) / (M 2 + — M 2 ) which measures the relative 
importance of the subleading terms in the expansion of the 7r + — ir° mass difference. 





Approx. 1 


Approx. 2 


Approx. 3 


z 

5 


0.899 
0.0177 


0.854 
0.0667 


0.852 
0.0683 



Table 1: Numerical results from the sum rule eq.(Q) corresponding to the central values 
of the physical parameters. Z and 5 are defined in the text. 

The contribution of the subleading terms is predicted to be positive and have a relative 
magnitude of 7%. This, of course, is in agreement with the upper bound that one obtains 
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Figure 6: Plot of the remainder part W^y(—s) in the sum rule integrand (see sec. 2) for 
the three successive approximations. 

from naive dimensional analysis of the low energy constants, which is 20%. What is the 
accuracy of this evaluation? We can identify three sources of error: l)the error coming 
from the uncertainties in the physical parameters that enter the calculation. 2) An error 
coming the chiral limit extrapolation and 3) an error associated with the assumption that 
A6 and As are constants, which is only an approximation. Concerning the first source 
of error, we have varied all the physical parameters independently and calculated the 
variation of the result for both approximations 2 and 3. The result is shown in table 2 
below. The parameters which are not shown like M p >, T p > induce very small errors. It is 
interesting that the individual errors are rather different in the two approximations. For 
instance, the error induced by F (here the chiral limit extrapolation error was included as 
well) or by Xq, As are significantly smaller in approx. 3 than in approx. 2. This does not 
imply that the third approximation has a smaller error, as it exhibits a greater sensitivity 
to the tail of the vector spectral function. If one simply adds all the errors one finds very 
closely the same number for the two approximations, respectively 7.3% and 7.5%. This 
is suggestive that both the central value and the error are just as reliably obtained from 
approx.2. In this approximation, we can also estimate the error due to the evaluation of 

o o 

the chiral limit values M p and T p . Varying the continuuum contribution in the set of sum 
rules as discussed in sec. 3.2 we obtain a small contribution of 0.2%. 
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parameters 


F 


A 6 


As 


M p 




P 


7 


M ai 


r ai 


error (2) 


0.9 


5.1 


1.0 


0.02 


0.04 


0.2 


0.05 






error(3) 


0.04 


3.4 


0.2 


0.01 


0.4 


1.4 


1.5 


0.2 


0.4 



Table 2: Percentage relative variation of the result forC/F 2 corresponding to the variation 
of the various input physical parameters within their error bars. The second and third 
lines of the table correspond to the calculation in approximation 2 and 3 respectively. 



The last uncertainty arises from the assumption made so far that A6 and As are con- 
stants which is only true at leading order in a s . This point can be investigated quantita- 
tively in the case of A6- Using the anomalous-dimension matrix provided in ref.(2*|], one 
can resum the leading logarithms and obtain, 



A 6 (s) 



647ra s (^) 
9 



1+ 9o^) lo s 



+(~o 6 a ( M ) - logo*)) 



| + 9o^) j_ 



4tt 



4tt 
-10/9 

+ 



9a g (/x) 
32vr 



-1/9 



l -fo^) + Ol(u) 



(47) 



In this expression, Oq and 0\ are the vacuum expectation values of the two operators 



A a - A a A a - A a 

O a & = {u^^—dd^^—u - u^—dd^—u) 

Oq = (Wy^dd^^u - u^y^dd^u) 



(48) 



with A a a color-space Gell-Mann matrix. In principle, in order to take the s-dependence 
correctly into account, one needs to know the values of both Oq and Oq. However, these 
two operators are not exactly on the same footing since Oq appears multiplied by one 
factor of a s more than Oq whenever the logarithm is not too large. Furthermore, Oq is 
suppressed in the large iV c limit. A plausible approximation, then, would be to ignore 
it altogether. Another plausible approximation is that of vacuum-saturation (2^], which 
yields the relation, 



(49) 



The energy dependence of Xq is very much suppressed in this approximation. As far as 
the integration over the circle in the complex plane is concerned (see eq. ([?])), we find that 
dropping the energy dependence is a very good approximation in any case, which does not 
generate an uncertainty in the determination of Xq(M t ) larger than 1%. One observes from 
eq.(|47|) that Xq is a steadily decreasing function of s which eventually goes to zero when s 
goes to infinity. Our construction can be seen as a procedure for smoothly matching the 



low-to- medium and the high energy regimes [16]. In this sense it is clear that one must 
choose Ag = Ag(so) where sq is the value of s where the asymptotic regime sets in, i.e. 
so must lie between M 2 and 2M 2 , say, as can be seen from Fig. 4. This determines the 
constant value of A6 to use within 2% approximately. Then, one must take into account 
the contribution of the logarithms in the high energy region of the sum rule integral. This 
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is found to introduce a rather small correction to the value of C which ranges from 0.4 to 
0.8% depending on the hypothesis made for 0\. In conclusion, we obtain that the overall 
relative error in the determination of the parameter C does not exceed 10%. 

Let us now consider the implication of this result for 0(p 4 ) low-energy parameters 
using the chiral expansion of the ir + — 7r° mass difference at this order 0, 

M^^-JL^^)) (5 0) 



167r 2 I fi 2 I F 2 \mci + m u 



with 



F k (v) = -2kl(ri + kl(v)+^) + Akl^) . (51) 

The 0(e 4 ) contribution which must technically be counted as 0(p i ) can be estimated to 
be numerically smaller by one order of magnitude than the 0(e 2 M 2 ) or the 0((m u — rrid) 2 ) 
ones and is neglected here. From this, and using m u /md = 0.55|4(]], one deduces, 

2.2 Z 7 + F k {M p ) = (-7.1 ±3.0)10~ 2 (52) 

which is our main result. For comparison, on the basis of naive dimensional analysis 
alone, one would obtain for the same quantity that it must lie in the range ±810~ 2 . The 
parameter l-j which appears in (|52|) is not very precisely known but a simple resonance- 



saturated sum rule gives an estimate |L0|] l-j ~ 0.710~ 2 . 
5. Conclusion 

To summarize, we have attempted an evaluation of the low-energy constant C with a 
controlled error, on the basis of the exact sum rule expression of Das et al.. The main prac- 
tical difficulty, which is present even if infinitely precise experimental data were available, 
lies in the necessity of extrapolating the integrand to the chiral limit. A calculational pro- 
cedure was proposed in which one first reconstructs the relevant current-current correlator 
in euclidian space making use of its smoothness properties together with the experimental 
determination of two asymptotic expansion parameters. An approximation scheme can 
be developped in which one includes spectral function components with higher and higher 
pion multiplicities. This expansion was argued to converge very rapidly such that, in 
practice, it is only necessary to include the one-pion and the two-pion components. The 
construction of the chiral limit makes use of recent work both on application of chiral per- 
turbation theory to the vector meson masses and of chiral calculations at two-loop order 
of current-current correlation functions. We have shown that under the assumption that 
the relative error on the asymptotic parameter A6 is of the order of 10% (this is the actual 
experimental error but it does not include the uncertainty stemming from the truncation 
of the OPE, which is more difficult to evaluate), one can determine the parameter C with 
an error of slightly less than 10% and deduce a meaningful estimate for a combination of 
subleading parameters k^. These parameters are primarily useful in calculations of radia- 
tive corrections at low energy. Another area where the computation of the photon loop 
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is of interest, is in relation with the K + — K° mass difference and the issue of Dashen's 
theorem violation. It is possible that the constraint obtained here may prove useful in this 
context as well. 
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